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a 

— 

prior  shape  parameter 

b 

- 

prior  scale  parameter 

b0 

- 
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- 
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- 
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Energy  Research  and  Development 

Administration 
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failures  per  hour 
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- 

prior  probability  density  function  of 

the  failure-rate  A 

f ( A | 0 failures) 

- 

posterior  probability  density  function  of 

the  failure-rate  A conditional  on  0 observed 

failures  in  a test  of  nt  unit-hours 

Y 

- 

posterior  assurance 

(1.0-y) 

- 

posterior  risk 

I (a,x) 

— 

incomplete  gamma  function 

jxya  ^e  ^dy 
0 

k 

- 

discrimination  ratio  (A^/A^) 

A 

- 

failure-rate 

xo 

- 

specified  failure-rate 

X1 

- 

criterion  failure-rate  (kAg) 

A*,  A* 

- 

generic  failure-rate 

iv 
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nt 
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?0 

P(-) 

P (x  failures | X ) 


P(0  failures) 


P ( 0 failures | X*<X<X*) 


P(x  failures) 


P(X#<X<X*|x  failures) 


50(l-pQ)th  percentile  of  a gamma  prior 

distribution;  lower  prior  limit 

Liquid  Metal  Fast  Breeder  Reactor 

Median  Failure  Rate 

Mean  Time  to  Failure 

number  of  test  units 

required  number  of  test  units 

unit-hours  of  testing 

required  unit-hours  of  testing 

prior  assurance 

probability 

conditional  probability  of  x failures  in 
a test  of  nt  unit-hours  with  an  under- 
lying failure-rate  of  X for  each  item 
unconditional  probability  of  passing  the 
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0 failures  in  a test  of  nt  unit-hours 
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failures  in  a test  of  nt  unit-hours 
unconditional  probability  of  x failures 
in  a test  of  nt  unit-hours 
conditional  probability  that  the  failure- 
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given  x failures  in  a test  of  nt  unit- 
hours 
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distribution 
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A BAYESIAN  ZERO-FAILURE  (BAZE)  RELIABILITY  DEMONSTRATION 


TESTING  PROCEDURE  FOR  COMPONENTS  OF  NUCLEAR 
REACTOR  SAFETY  SYSTEMS 

Abstract 

A Bayesian-Zero-Failure  (BAZE)  reliability  demonstration 
testing  procedure  is  presented.  The  method  is  developed  for  an 
exponential  failure-time  model  and  a gamma  prior  distribution 
on  the  failure-rate.  A simple  graphical  approach  using  percentiles 
is  used  to  fit  the  prior  distribution.  The  procedure  is  given 
in  an  easily  applied  step-by-step  form  which  does  not  require 
the  use  of  a computer  for  its  implementation.  The  BAZE  approach 
is  used  to  obtain  sample  test  plans  for  selected  components  of 
nuclear  reactor  safety  systems. 


II 


vii 


I.  INTRODUCTION 

Most  government  and  military  contracts  for  hardware  develop- 
ment include  a numerical  reliability  requirement  in  the  specifications. 
For  example,  a certain  recent  contract  required  an  overall  system 
failure-rate  no  larger  than  11.64  x 10  ^ f/h.  In  addition,  most 
contracts  require  quantitative  demonstrated  assurance  that  such  a 
requirement  has  been  met.  MIL-STD-781C+  provides  a standard  which 
can  be  used  to  demonstrate  such  a requirement  for  times-to-f ailure 
that  are  exponentially  distributed.  The  standard  may  be  used  for 
preproduction  (qualification)  tests,  as  well  as  production 
reliability  acceptance  (demonstration)  tests.  A typical  reliability 
demonstration  statement  is  that  a failure-rate  requirement  of  X 
failures/h  be  demonstrated  with  Y%  confidence. 

The  purpose  of  this  report  is  to  develop  a Bayesian 
reliability  demonstration  testing  procedure  for  exponentially 
distributed  failure  times  which  can  be  easily  and  effectively 
used  to  demonstrate  component/subsystem/system  reliability 
conformance  to  stated  requirements.  The  procedure  will  also 
be  used  to  develop  suggested  test  plans  for  various  components 
used  in  nuclear  power  reactor  safety  systems.  This  procedure 
may  be  used  to  verify  specified  and  projected  component 
failure-rates  in  LMFBR  safety  systems.  However,  the  procedure 
is  a general  one  and  its  use  is  not  restricted  to  nuclear  power 

+ MIL- STD- 7 8 1C  "Reliability  Qualification  and  Production 
Acceptance  Tests:  Exponential  Distribution,"  Washington,  D.C.: 

U.S.  Government  Printing  Office,  (Draft),  August,  1976. 
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safety  systems.  It  may  also  be  used  to  demonstrate  reliability 
for  such  equipment  categories  as  ground  equipment,  shipboard 
equipment,  avionic  equipment,  weapons  systems,  and  surveillance 
equipment. 

Over  the  past  two  decades  numerous  classical  reliability  de- 
monstration testing  methods  have  been  devised  for  various  failure- 
time distributions.  Classical  demonstration  test  plans  for  com- 
ponents having  an  exponential  failure-time  distribution  (constant 
failure-rate)  may  be  fixed  time  (Type  I Censoring) , fixed  number 
of  failures  (Type  II  Censoring) , or  sequential  tests.  In  addition, 
such  tests  may  be  conducted  either  with  or  without  the  replacement 
of  failed  items  when  failures  occur  during  the  test.  One  example 
of  a classical  procedure  is  MIL-STD-781C  which  gives  various  test 
procedures  derived  under  the  assumption  of  a constant  failure- 
rate.  These  tests  are  based  on  various  levels  of  producer's  and 
consumer's  risk,  design  ratios,  and  the  confidence  level  of  the 
test.  In  these  tests,  a simple  statistical  hypothesis  of  the 
form 

H:  A=Aq  (specified  failure-rate) 

A:  A=A^  (maximum  acceptable  failure-rate) 

is  tested,  where  A is  the  (unknown)  failure-rate  of  the  device  and 
Ai>Aq-  The  producer's  risk  is  the  probability  that,  if  H is 
true,  A will  be  accepted,  while  the  consumer's  risk  is  the 
probability  that,  if  A is  true,  H will  be  accepted.  The  design 
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The  text  by  Mann 


ratio  is  defined  as  the  ratio  of  A^  to  Ag. 

Schafer,  and  Singpurwalla  (1974)  gives  an  excellent  discussion 
of  the  basis  upon  which  these  and  other  classical  test  plans 
have  been  developed. 

In  practice,  it  is  often  the  case  that  the  reliability 
parameters  of  interest,  such  as  MTTF , reliable  life,  failure- 

O 

rate,  etc. , most  realistically  should  be  treated  as  random 
variables  and  not  as  constant  values.  The  statistical  distri- 
bution which  expresses  the  true  underlying  variation  in  the 
parameter  is  called  the  prior  distribution  of  the  parameter 
(when  treated  as  a random  variable) . This  approach  has  been 
taken  in  previous  reactor  safety  analyses,  such  as  the  Rasmussen 
study.*  In  that  study,  the  prior  distribution  for  the  reactor 
component  failure-rates  was  taken  to  be  the  log-normal  distri- 
bution.** Such  an  approach  permits  the  reliability  quantity  of 
interest  to  vary  randomly  due  to  such  factors  as  environmental 
effects,  plant-to-plant  differences,  maintenance  effects,  and 
different  operational  demands. 

By  treating  the  parameter  as  a random  variable,  a Bayesian 
approach  may  be  considered.  The  main  advantage  of  the  Bayesian 
approach  is  that  the  resulting  estimates  are  computed  from  all 
available  information,  and  not  just  narrowly  defined  test  data 

*WASH  1400  Appendix  III 

**Ibid. , p.  11-40 
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of  precise  content.  Rather,  there  exist  two  sources  of  information 
regarding  the  Bayesian  procedure.  One  source  of  information,  the 
so-called  "prior  information,"  expresses  the  sum  total  of  engineer- 
ing judgment  and  belief  concerning  the  underlying  prior  distribu- 
tion of  the  parameter  of  interest.  It  is  precisely  this  distri- 
bution which  expresses  the  innerent  variability  of  the  parameter 
itself.  The  other  source  of  information  is  the  statistical  model 
used  to  describe  either  the  time-to-f ailure  data  or  the  test  re- 
sults themselves.  Both  sources  of  information  are  combined  via 
Bayes  Theorem  [see  Waller  and  Martz  (1975) ] to  produce  a statement 
sucn  that  the  probability  that  the  failure-rate  does  not  exceed 
the  specified  value  is  Y%. 

Such  an  approach  is  particularly  applicable  for  deriving  test 
plans  for  demonstrating  the  component  failure  rates  of  proposed 
nuclear  reactors.  The  reason  for  this  is  that  failure  data  are 
becoming  available  for  similar  components  in  use  in  existing 
power  reactor  systems  througnout  the  world.  These  data,  which  are 
continuously  being  compiled  and  reported  in  numerous  data  base 
systems,  represent  the  "prior  information"  for  similar  components 
to  be  used  in  advanced  reactor  systems,  such  as  LMFBR  systems. 

The  resulting  Bayesian  test  plans  are  generally  resource- 
effective,  due  to  the  use  of  all  available  information  and 
judgment  concerning  the  parameter  of  interest.  If  the  prior 
information  supports  an  adequately  reliable  component,  then  less 
testing  will  usually  be  required  compared  to  the  classical  case. 
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If  the  opposite  is  true,  then  more  testing  may  be  required.  A 
general  introduction  to  the  use  of  Bayesian  methods  in  reliability 
is  given  by  Waller  and  Martz  (1975),  (19'<  Oa)  , (1976b). 

A brief  review  of  Bayesian  reliability  demonstration  proce- 
dures is  given  in  Section  II.  A procedure  for  choosing  the  prior 
distribution  is  presented  in  Section  III.  The  BAZE  procedure  is 
developed  in  Section  IV.  Section  V contains  an  example  illustra- 
tion of  the  method,  as  well  as  an  examination  of  the  sensitivity 
of  the  results  to  the  chosen  prior  distribution.  Appropriate 
prior  distributions  are  fitted  to  various  components  used  in 
nuclear  power  systems  in  Section  VI.  Section  VII  presents  a 
selection  of  sample  BAZE  test  plans  for  these  reactor  components. 

II.  BAYESIAN  RELIABILITY  DEMONSTRATION  TESTING 

The  state  of  the  art  of  Bayesian  reliability  demonstration 
test  procedures  will  now  be  reviewed.  One  of  the  earliest  re- 
ferences to  Bayesian  reliability  demonstration  plans  is  that  of 
Bonis  (1966) . Since  that  time  numerous  Bayesian  schemes  have 
been  developed.  Easterling  (1970)  presented  a somewhat  modified 
Bayesian  demonstration  procedure.  Schafer  and  Singpurwalla  (1970) 
developed  a sequential  Bayes  procedure  for  obtaining  required 
test  plans.  Schafer  (1969)  , (1971) , and  (1973)  has  considered 

three  types  of  Bayesian  plans:  (1)  Bayesian  fixed  time  tests, 

(2)  mixed  Bayesian/classical,  and  (3)  Bayesian  sequential 
tests.  Following  along  these  same  lines,  Goel  et  al . , (dates 


unknown)  developed  Bayesian  plans  for  slightly  different  criteria. 
Blumenthal  (1973)  has  also  developed  Bayesian  test  plan  procedures. 
Guild  (1968),  (1973)  has  developed  what  he  refers  to  as  "median 
failure  rate"  (MFR)  reliability  demonstration  plans.  Other 
Bayesian  plans  have  also  been  considered  by  Balaban  (1969) , (1975) 

and  Ramos  (1970) . Joglekar  (1975)  discusses  several  of  these 
Bayesian  testing  schemes.  Recently,  Goel  and  Joglekar  [1976] 
have  prepared  a comprehensive  account  of  the  state  of  the  art 
of  Bayesian  reliability  acceptance  sampling.  This  five-part 
series  provides  an  excellent  introduction  to  the  subject. 

One  of  the  major  problems  with  most  Bayesian  test  plans  is 
the  relative  difficulty  in  obtaining  a desired  plan  in  practice. 
This  is  due  to  the  presence  of  additional  prior  parameters,  as 
well  as  the  relative  complexity  of  the  method.  For  example,  most 
of  the  above  Bayesian  plans  are  derived  for  a constant  failure- 
rate  model  and  a gamma  prior  distribution  on  the  failure-rate 
[see  Waller  and  Martz  (1975)]. 

The  gamma  prior  distribution  is  the  natural  conjugate  prior 
distribution  for  the  constant  failure-rate  model.  Schafer  (1969) 
investigated  data  from  32  different  equipments  and  found  that  in 
29  cases  a gamma  prior  distribution  adequately  fit  the  data. 

Others  have  likewise  observed  the  suitability  and  versatility 
of  the  gamma  prior  distribution.  For  these  reasons  it  is  also 
considered  here.  However,  it  is  frequently  not  an  easy  task  to 
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identify  an  appropriate  gamma  prior  distribution  and,  once  this 
has  been  done,  to  obtain  the  required  test  plan.  The  choice  of 
test  plan  criteria,  e.g.,  the  consumer  and  producer  risks  that 
are  to  be  controlled,  can  also  complicate  the  determination  of 
the  required  plan.  ^rtain  criteria  yield  plans  that  are  simple 
to  determine.  One  of  the  easiest  of  the  above  Bayesian  proce- 
dures to  use  in  practice  is  that  given  by  Guild  (1973) . A more 
general  version  of  this  procedure  is  developed  in  Section  IV 
for  use  here.  An  earlier  version  of  this  BAZE  procedure  was 
presented  by  Martz  and  Waller  (1976c) . The  current  BAZE  procedure 
contains  a more  useful  and  practical  procedure  for  selecting  the 
prior  distribution,  which  has  been  incorporated  into  the  BAZE 
procedure  itself.  The  new  procedure  also  contains  a simple 
method  for  examining  the  sensitivity  of  the  resultant  BAZE  test 
plan  to  the  chosen  prior  distribution.  This  serves  to  make  the 
method  more  useful  in  practice. 

Three  somewhat  distinctive  aspects  of  the  procedure  should 
be  mentioned  before  we  begin  the  development.  First,  the  criterion 
upon  which  the  procedure  is  based  is  simple,  pertinent,  and  easy 
to  grasp.  Second,  the  fitting  of  the  prior  distribution  is  an 
integral  part  of  the  procedure  and  is  based  on  the  use  of  informa- 
tion regarding  two  percentiles  of  this  distribution.  Third,  the 
procedure  is  straightforward  and  easy  to  apply  in  practice, 
with  only  a few  simple  graphs  and  tables  and  pocket  calculator 
required.  Together  these  provide  a useful  Bayesian 
procedure  for  a large  variety  of  applications. 
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III.  SELECTING  A GAMMA  PRIOR  DISTRIBUTION 


We  assume  that  the  time-to-f ailure  of  interest  is  an 
exponentially  distributed  random  variable  with  failure-rate 
parameter  A.  For  Bayesian  analyses  in  this  model,  the  family 
of  gamma  distributions  with  probability  density  functions 
given  by 

f(X)  = rTIT  Xa_1e“bX,  X,  a,  b > 0 (1) 

provides  conjugate  prior  models  for  X.  In  practice,  an  engir^er 
must  select  a member  of  this  family  as  the  prior  distribution 
to  be  used  in  determining  the  BAZE  test  plan  to  be  discussed  in 
the  next  section.  The  selection  of  a particular  prior  distribu- 
tion is  accomplished  by  identifying  values  for  the  prior  shape 
parameter  a and  prior  scale  parameter  b.  The  parameter  a can 
be  further  interpreted  as  the  number  of  pseudo  failures  in  a 

prior  life  test  of  duration  b pseudo  hours.  The  mean  and 

2 

variance  of  A are  given  by  (a/b)  and  (a/b  ) , respectively. 

Some  additional  benefits  in  using  a gamma  prior  distribution 
are  as  follows: 

a.  The  two  parameters  give  sufficient  flexibility  to 
model  a variety  of  shapes  of  prior  distributions  likely  to 
be  encountered  in  practice.  The  following  indicate  the  possible 
shape  characteristics  of  a conjugate  gamma  prior  distribution. 
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b>  0 


L-shaped 

Unimodal 

or 

Exponential 

with 

mode  at 

decreasing 

b(a-l) 

0<a<l 

a=l 

a>l 

b.  The  positive  skewness  can  account  for  general  behaviors 
of  assessed  data  in  which  less  likely  but  large  deviations  may 
occur  (such  as  abnormally  high  failure-rates  due  to  batch  defects, 
environmental  degradation,  and  other  outlier  causing  effects) . 

c.  In  practice,  the  gamma  family  often  satisfactorily  fits 
observed  data  [see  Schafer  (1969)]. 

d.  The  positively  skewed  nature  of  the  gamma  family  provides 
a protective,  positive-type  bias  which  is  retained  when  the 
distribution  is  propagated  by  means  of  a Bayesian  analysis. 

A simple  method  for  determining  values  for  a and  b will  now 
be  described.  The  method  requires  an  engineer  to  provide  upper 
and  lower  percentile  values.  Once  these  are  given,  a simple 
graphical  or  table  look  up,  in  addition  to  a few  simple  calcula- 
tions, yields  the  corresponding  values  of  a and  b. 

The  engineer  must  provide  two  values  of  the  failure-rate  X, 
referred  to  as  the  lower  prior  limit  (LL)  and  upper  prior  limit 
(UL) , such  that 
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P ( A < LL) 


P ( A > UL) 


(1.0-p0)/2 


(2) 


where  Pq  is  required  to  be  equal  to  one  of  the  values  0.95,  0.90, 
or  0.80,  and  where  LL  < UL.  That  is,  LL  and  UL  are  specified 
such  that  there  is  an  equal  50(1.0-pQ)%  chance  that  the  true 
(unknown)  failure-rate  is  either  less  than  LL  or  greater  than  UL, 
respectively.  Thus  LL  and  UL  are  the  respective  50(1.0-PQ)th 
and  50(1.0+PQ)th  percentiles  of  the  prior  gamma  distribution. 

For  example,  suppose  that  an  engineer's  best  prior  judgement  or 
belief  is  that  P(A  < 1.0x10  ^f/ii)-b%  and  that  P(A  > 1.0x10  ^f/h)=5%. 
Thus  (1.0-Pg)/2  = 0.05,  Pq  = 0.90,  LL  = 1.0x10  ^f/h,  and  UL  = 

1.0x10  ^f/h.  The  quantity  100pg%  is  the  prior  assurance  that 
the  interval  (LL,  UL)  contains  the  failure-rate  of  interest. 

Since  engineers  are  increasingly  becoming  accustomed  to  working 
with  5%  error  probabilities,  it  is  likely  that  Pq  = 0.90  will 
normally  be  used.  However,  80%  and  95%  prior  assurances  can 
also  be  used.  For  example,  in  the  Rasmussan  study  (WASH-1400), 

90%  prior  assurance  was  considered.  If  the  prior  assurance  is 
free  to  be  selected,  it  is  recommended  that  Pq  = 0.90  be  used. 

In  this  case,  LL  and  UL  become  the  lower  and  upper  prior  5% 
bounds,  respectively,  on  the  failure-rate. 

A mathematical  justification  of  the  procedure  to  be  described 
is  given  in  Appendix  B.  A step-by-step  outline  of  the  procedure 
is  as  follows: 
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Step  1:  Specify  the  values  of  LL,  UL,  and  Pg  = 

0.90,  or  0.95  that  represent  the  totality 
of  your  best  judgement  and  belief  about  the 
failure-rate  A of  interest.  These  values  are 
selected  in  accordance  with  (2). 

Step  2:  Compute  the  value  of  log^g  (UL/LL) . 

Step  3:  For  the  value  of  Pg  chosen  in  Step  1 and  the 

value  of  log1Q  (UL/LL)  calculated  in  Step  2, 
read  the  required  value  of  shape  parameter  a 
from  Figure  A1  (see  Appendix  A) . 

Step  4 : For  the  value  of  Pg  from  Step  1 and  for  the  value 

of  a found  in  Step  3,  read  the  value  of  bg  from 
Figure  A2  (see  Appendix  A) . 

Note:  Table  A2  (see  Appendix  A)  may  be  used  in 

lieu  of  Figure  A2  to  obtain  bg,  depending  upon 
which  is  more  convenient  to  use.  If  necessary, 
interpolate  in  Table  A2 . 

Step  5:  For  the  value  of  LL  from  Step  1 and  the  value  of 

bg  from  Step  4,  calculate  the  required  value  of 
the  scale  parameter  b according  to 

b = bg (1 . 0xl0-6f/h) /LL. 

Let  us  illustrate  this  procedure  by  means  of 
an  example. 
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Example:  For  a certain  component  of  interest,  suppose  it  is 

is  believed  that  the  failure-rate  A is  such  that 
P ( A < 1 . 0xl0~7f/h)  = 5%  and  P(A  > 1.0xl0-5f/h)  = 5%. 

It  is  required  to  identify  the  particular  gamma  distribution 
which  is  consistent  with  this  belief. 

Step  1:  LL  = 1.0x10  7f/h,  UL  = 1.0x10  ^f/h,  and  p^  = 0.90. 

Step  2:  log1Q (UL/LL)  = log10(102)  = 2.0. 

Step  3:  From  Figure  Al,  for  pQ  = 0.90  and  log^Q (UL/LL)  = 

2.0,  we  find  a = 0.84. 

Step  4:  For  a = 0.84  and  pQ  = 0.90,  Table  A2  yields 

bQ  = 2.6723xl04h. 

Step  5:  The  required  scale  parameter  b becomes 

b _ (2.6723xl04h) ( 1 . 0xl0~6f /h) 

( 1 . 0xl0~7f/h) 

= 2 . 6723xl05  h. 


was  done  to  insure  that  the  lower  tail  area  will  always  be 
( l~Pg) /2 , while  the  upper  tail  area  may  depart  somewhat  from 
the  desired  value  (l-pQ)/2.  We  chose  to  hold  the  lower  tail 
area  fixed  because  of  the  positively  skewed  nature  of  the  gamma 
distribution . 

One  final  note  concerns  the  usefulness  of  Figures  A1  and 
A2  in  practice.  The  effective  range  of  values  of  a considered 
in  Figures  Al  and  A2  is  between  0.25  and  10.0.  Experience  with 
fitting  gamma  prior  distributions  to  failure-rate  data  indicates 
that  this  range  should  contain  nearly  all  situations  likely  to 
be  encountered  in  practice.  This  range  is  consonant  with  ratios 
of  UL  to  LL  roughly  between  0.3  and  4.0  orders  of  magnitude 
(powers  of  10) . 


IV.  BAYESIAN  ZERO-FAILURE  (BAZE)  RELIABILITY  DEMONSTRATION 
TESTING 


The  BAZE  reliability  demonstration  procedure  was  developed 
by  Martz  and  Waller  (1976c) . This  procedure  is  given  here  in  an 
expanded  form  which  includes  the  procedure  for  fitting  a gamma 
prior  distribution,  described  in  the  preceding  section,  and 
other  features  as  well. 

This  section  describes  how  to  construct  and  apply  Bayesian 
fixed  time  demonstration  test  plans  of  the  replacement  type, 
called  BAZE  plans,  for  systems/subsystems/components  having  a 
constant  failure-rate.  The  BAZE  procedure  is  appropriate  for 
testing  time-dependent  chance  failure  mechanisms. 
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To  begin,  consider  a device,  henceforth  referred  to  as  a 
"component,"  having  an  exponential  failure  time  distribution 
with  failure-rate  A.  Thus,  the  failure-time  random  variable  X 
of  this  component  is  assumed  to  follow  the  well-known  exponential 
probability  density  function  given  by 

f (x| A)  = Ae~Xx,  x > 0,  A > 0.  (3) 

As  mentioned  earlier,  it  is  assumed  here  that  the  prior  distribu- 
tion of  A is  the  natural  conjugate  gamma  distribution  with 
probability  density  function  given  by 

f(A)  = rTil  xa”le"bX'  x»  a'  b > 0,  (4) 

where  a is  the  prior  shape  parameter  and  b is  the  prior  scale 
parameter. 

The  test  plans  considered  here  assume  that  n identical 
components  are  tested  each  for  a prespecified  length  of  time  t, 
the  test  duration.  The  quantities  n and  t are  to  be  determined 
consistent  with  the  following  statistical  test  criterion.  The 
test  criterion  is  as  follows:  if  no  failures  occur,  the  test 

is  passed,  while  if  one  or  more  failures  occur,  the  test  is 
failed.  Thus  the  test  is  terminated  either  at  the  prespecified 
test  time  t or  at  the  time  of  the  first  component  failure,  which- 
ever occurs  first.  Such  "zero-failure"  test  plans  usually 
require  the  smallest  unit-hour  test  combination  nt  for  a stated 
confidence,  and  are  thus  test  resource-effective.  In  addition. 
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by  restricting  consideration  to  zero-failures,  such  test  plans 
are  easy  to  obtain.  Now  the  probability  of  obtaining  exactly 
zero  failures  during  the  test  is  given  by 

P (0  failures|A)  = e nt\  (5) 

The  posterior  distribution  of  A is  also  a gamma  distribution 
with  scale  parameter  (b  + nt)  and  shape  parameter  a.  Thus, 
conditional  on  zero  failures  in  nt  unit-hours  of  testing,  the 
posterior  probability  density  function  of  A becomes 

£ ( A | 0 failures)  - »a'V  lb+nt)  X . (6) 

In  order  to  find  the  required  unit-hour  test  combination 
nt,  some  criterion  regarding  the  desired  confidence  level  of 
the  demonstration  test  must  be  given.  The  plans  presented  here 
satisfy  the  posterior  risk  criterion  given  by 

P ( A < kig  |0  failures)  = y,  (7) 

where  P(-)  is  a probability  function.  Here  Ag  is  the  specified 
failure-rate . If  we  define  A^  = kAg,  t*ien  xi  referred  to  as 
the  test  criterion  failure-rate  and  k is  known  as  the  discrimination 
ratio.  The  test  criterion  in  (7)  is  interpreted  as  follows.  In 
a test  that  is  passed,  i.e.,  zero  failures  occur,  the  probability 
is  100y%  that  the  component  failure-rate  does  not  exceed  (kAg) . 
Here(l.O-y)  will  be  referred  to  as  the  posterior  risk  and  y will 
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be  referred  to  as  the  posterior  assurance.  The  word  "posterior 


denotes  that  the  assurance  pertains  to  tests  which  have  been 
passed.  Of  course,  values  of  y and  k are  required  in  order  to 
determine  the  required  test  plan.  More  will  be  said  about  the 
selection  of  these  values  later. 

Recall  that  the  test  procedure  requires  that  n items  be 
placed  on  life  test  for  t hours.  If  no  failures  occur,  the  test 
is  passed.  Now,  if  the  test  is  passed,  it  may  be  claimed  that  a 
failure-rate  not  exceeding  (kXQ)  has  been  demonstrated  with 
100y%  posterior  assurance.  If  a single  failure  occurs,  the 
test  is  failed,  and  the  forgoing  claim  cannot  be  made. 

The  BAZE  procedure  described  requires  specification  of 
values  for  the  following  five  quantities: 

— LL( lower  prior  limit) 

— UL (upper  prior  limit) 

- Pg(the  prior  assurance) 

— X^(the  criterion  failure-rate) 

— y (the  posterior  assurance). 

It  is  noted  that  the  criterion  failure-rate  X^  may  be  equal 
to  the  specified  failure-rate  Xq.  In  this  case  k=l;  otherwise, 
k^l.  The  procedure  is  developed  by  writing 

,,  ^ -i  > r*  (b+nt)a  ,a*l  -(b+nt)Xj, 

p ( X < X * I 0 failures)  = J - — . — X e dX 

0 1 ia; 

. I (a , [b+nt] X * ) . 

r(a)  ' 1 
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where  I(a,x)  is  the  widely  studied  incomplete  gamma  function 
defined  by 


I(a,  x) 


rx  a-1  -y, 

J Y e Jdy . 

0 


Tables  and  computer  routines  for  evaluating  this  function  are 
widely  available  for  use  in  our  development.  Hence,  when  LL , 
UL,Pq,Y/  and  ^ are  specified,  the  step-by-step  procedure  for 
obtaining  the  required  BAZE  test  plan  is  as  follows: 


Step  1: 


Step  2: 


Step  3: 


Step  4 : 


For  the  specified  values  of  LL  and  UL,  compute 
the  value  of  log^g (UL/LL) . 

For  Pq  = 0.80,  0.90,  or  0.95,  and  the  value  of 
logiO (UL/LL)  from  Step  1,  obtain  the  value  of 
the  prior  shape  parameter  a from  Figure  A1 
(Appendix  A) . 

For  Pg  and  the  value  of  a from  Step  2,  obtain 
the  value  of  bg  from  either  Figure  A2  or  Table  A2 . 
For  the  value  of  LL  and  bg  from  Step  3,  calculate 
the  value  of  the  prior  scale  parameter  b,  in 
appropriate  time  units,  according  to 


b = bg  ( 1 . 0xl0~6f/h) /LL. 

Step  5:  Obtain  the  value  of  0^  from  Table  Al  (Appendix  A) 

for  the  value  of  a found  in  Step  2 and  y.  Note : 
Table  Al  may  be  used  directly  for  a = 0.0001 


(0.0001)  0.01  (0.001)  0.10  (0.01)  1.0  (0.1)  5.0 


(0.5)  10.0  (1.0)  50.0  and  y = 0.99,  .975,  .95, 

.90,  .85,  .80,  .75,  .70,  .60,  .50.  For  other 
values  of  a and/or  y,  either  interpolate  in  Table 
A1  or  solve  the  equation  given  by 

/ Y Xa-1  e~A  dX  - y r (a)  = 0 
0 

for  0 . It  is  mentioned  here  that  in  constructing 
Table  A1  the  incomplete  gamma  function  in  the  above 
equation  was  numerically  calculated  by  use  of  the 
code  INCGAM,  written  by  D.E.  Amos  and  S.L.  Daniel 
of  Sandia  Laboratories,  Albuquerque,  NM,  November 
1974.  The  above  equation  was  solved  on  a CDC  6600 
computer  by  use  of  the  root-solving  code  ZEROIN, 
written  by  L.F.  Shampine  and  H.A.  Watts,  also  of 
Sandia  Laboratories,  September  1970. 

Step  6:  With  X^  and  the  6^  value  from  Step  5,  and  the  value 

of  b from  Step  4,  solve  for  the  required  BAZE  unit- 
hours  of  test  (nt)Q  given  by 

(nt)  Q = (9  - bX^/X^ 

Note:  Negative  values  of  (nt)Q  can  occur.  A nega- 

tive value  of  (nt)g  can  be  interpreted  as  a demon- 
stration of  the  failure-rate  x^,  at  the  stated 
posterior  assurance  level,  without  the  need  for 
further  testing.  This  situation  occurs  whenever 


the  prior  distribution  satisfactorily  meets  the 
assurance  requirement. 

! 

Step  7:  Calculate  the  sensitivity  of  the  final  plan  and 

the  posterior  assurance  to  errors  in  the  prior 
assurance,  as  well  as  the  sensitivity  of  the 
final  plan  to  errors  in  the  posterior  assurance 
and  criterion  failure-rate  (cf  Section  V.  Example) . 

Step  8:  Identify  the  final  plan  to  be  used  based  on  the 

results  of  the  sensitivity  analysis  conducted  in 
Step  7 and  the  unconditional  probability  of  passing 
the  test  [from  (10)  or  (11)  below]. 

Step  9:  The  required  test  duration  t^  and  number  of  test 

units  n^  is  given  by  any  pair  of  values  satisfying 
n^tg  = (ntj^,  where  (nt)g  is  the  required  unit-hours 
of  test  from  Step  8 and  nQ  is  a positive  integer. 

The  values  ng  and  t^  are  selected  by  outside  con- 
siderations, such  as  test  time  constraints  and  the 
number  of  test  units  available. 

It  is  noted  here  that  if  a = 1,  b = 0,  and  y = 0.50,  then 
the  BAZE  test  plan  is  exactly  the  same  as  the  classical  test  plan. 
If  a = 1 and  y = 0.50,  then  the  BAZE  test  plan  will  always  require 
less  unit-hours  of  testing,  depending  upon  the  magnitude  of  b. 

- For  a specified  fixed  value  of  a,  the  required  BAZE  unit-hours 

of  test  decreases  as  y decreases,  b increases,  or  k increases. 
Consequently,  for  given  values  of  a and  b,  an  opportunity  is 
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present  for  the  test  designer  to  trade  between  decreasing  testing 
costs  [decreasing  (nt)^]  and  decreasing  test  assurance  [decreasing 
y and/or  increasing  A^].  Such  tradeoffs  are  illustrated  in 
Section  V. 

Suppose  that  high-reliability  components  with  failure-rate 
on  the  order  of  magnitude  of  10  ^ f/h  are  being  considered.  If 
the  prior  mean  is  in  this  range,  then  it  is  true  that  as  the 
spread  between  LL  and  UL  increases  both  a and  b will  generally 
decrease.  This  is  seen  in  Figures  A1  and  A2 . In  fact,  in 
situations  such  as  this,  a will  frequently  be  less  than  one. 

This  situation  occurs  whenever  the  prior  variance  is  quite  large, 
i.e.,  whenever  the  prior  distribution  is  diffuse  [see  Waller  and 
Martz  (1975)].  Such  situations  frequently  occur  in  reliability 
and  this  fact  has  motivated  the  fine  grid  of  a values  less  th_.i 
one  considered  in  Table  Al . 

A quantity  of  particular  interest  to  the  producer  is  the 
unconditional  probability  of  passing  the  test  when  using  a test 
plan  with  (nt)  unit-hours  of  test.  The  unconditional  probability 
of  not  passing  the  test  is,  in  some  sense,  the  "producer's  risk" 
of  the  BAZE  procedure.  The  probability  of  passing  the  test  must 
be  sufficiently  large  in  order  that  the  producer  be  willing  to 
conduct  the  test.  This  probability  also  conveys  to  the  consumer 
the  likelihood  that  the  required  posterior  assurance  will  be 
realized.  This  probability  is  given  by 
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P(Passing  the  Test)  e P(0  failures)  = [b/tb+nt)]a.  (10) 

Related  to  this  quantity  is  the  conditional  probability  of 
passing  the  test  when  it  is  known  that  the  true  (unknown)  failure- 
rate  A lies  within  a given  interval  [A*,  A*],  where  0<A*<A*<<». 

In  this  case  we  have 


P(Passing  the  Test | A *< A < A * ) = 


(_b_T 

I (a, [nt+b] A*) -I (a, [nt+b] A*) 

l b+nt  J 

I (a , b A * ) -I (a, bA*) 

t 

(ID 


where  I(a,x)  is  defined  in  (9).  It  is  also  noted  that,  if  A*  = 0 
and  A*  = °°,  then  the  conditional  probability  of  passing  the  test 
given  in  (11)  reduces  to  the  unconditional  probability  given  in 
(10).  In  practice,  an  interval  [A*,  A*]  which  is  certain  to  contain 
the  failure-rate  can  frequently  be  identified.  If  this  can  be  done, 
then  (11)  should  be  used  in  place  of  (10) . 

What  posterior  assurance  do  we  have  about  the  failure-rate 
if  one  or  more  failures  occur  during  the  test,  i.e.,  if  the  test 
is  failed?  Suppose  that  failed  items  are  replaced  as  they  occur 
during  the  test.  Then 


P (x  failures | A ) 


-nt  A x 

e (nt  A ) 


, x—  0,  1,  ..., 


(12) 


and  the  posterior  probability  density  function  of  A becomes 

a+x 


f ( A x failures) 


(b+nt) ' 
r (a+x) 


a+x-1  - (b+nt) A 

A G 


(13) 


Now,  for  any  specified  interval  [A*,  A*],  where  0<A*<A*<°°,  we  have 
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P(A*<A<A*|x  failures)  = 

r (a+x)  £ I^a  + x'  tb  + nfc]A*)  ~ Ka  + x*  [b  + nt]A*)  j. 

In  particular,  we  have 


(14) 


P(A<kA  |x  failures)  = I (a  + x,  [b  + nt]kA0)/I’(a  + x)  . (15) 

It  is  observed  that,  if  x = 0,  then  (15)  is  equal  to  the  specified 
posterior  assurance  y.  For  the  criterion  failure-rate  ( k A ^ ) , as 
x increases  (15)  becomes  smaller  than  y.  Thus,  as  more  failures 
occur,  we  have  less  posterior  assurance  about  the  failure-rate  not 
exceeding  the  criterion  value. 

Also,  the  unconditional  probability  of  obtaining  exactly  x 
failures  in  a test  of  nt  unit-hours  duration  is 


P(x  failures)  = 


ba (nt) xr (a  + x) 


r (a)  r (x  + 1)  (nt  + b) 


a + x 


(16) 


The  statistical  performance  characteristics  of  the  chosen 
plan  are  completely  summarized  by  means  of  the  posterior  operating 
characteristic  (POC)  curve.  This  curve  is  obtained  by  plotting 
P ( A < A * | 0 failures)  as  a function  of  A*.  Unlike  classical  OC  curves, 
the  POC  curve  is  a cumulative  distribution  function.  This 
probability  may  be  computed  from  (8) . 

We  cannot  emphasize  enough  that  both  the  consumer  and 
producer  must  be  willing  to  pay  the  price  for  increasing  assurance 
of  small  failure-rates  by  increasing  the  unit-hours  of  testing. 

This  will  be  illustrated  in  the  example  in  the  next  section. 
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V.  EXAMPLE 

In  order  to  fully  illustrate  the  BAZE  procedure,  consider 
the  following  example.  Consider  a certain  component  whose 
random  failure-rate  is  required  to  be  demonstrated.  How  many 
unit-hours  of  testing  with  no  countable  failures  are  required 
in  order  to  be  able  to  claim  that  P (A<2 . 0xl0-8f /h)  > 0.70, 
after  the  test  has  been  passed?  From  past  experience  and  engi- 
neering judgement,  suppose  it  is  believed  that 

P{X<8.5xlO-8f/h}  = P{ A>4 .8xi0-6f/h}  = 5%. 

Thus  X1  = \Q  = 2 . 0xl0_6f/h,  y = 0.70,  LL  = 8.5xl0~8f/h,  UL  = 

4.8x10  8f/h,  and  p^  = 0.90.  Following  the  step-by-step  proce- 
dure in  the  preceding  section  yields  the  following  results: 

Step  1:  log1Q (UL/LL)  = 1.75. 

Step  2:  For  log^Q (UL/LL)  = 1.75  and  p^  = 0.90,  we  read 

a = 1.0  from  Figure  A1 . 

Step  3:  For  p^  = 0.90  and  a = 1.0,  we  obtain  b^  = 5.1293xl04h 

from  Table  A2 . 

Step  4:  For  LL  = 8.5x10  8f/h  and  b^  = 5.1293xl04h,  we 

calculate 

b = (5.1293xl04h)  (1. 0xl0~6f/h)/( 8.5x1 0~8f/h) 

= 603 . 44xl03h. 

Step  5:  For  a = 1.0  and  y = 0.70,  we  obtain  0^  ^ = 1.203973 

from  Table  Al. 
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Step  6:  For  A ^ = 2.0x10  ^f/h,  6g  = 1.203973,  and 

b = 0.60xl0^h,  the  required  BAZE  unit-hours 
of  test  (nt)g  are  calculated  to  be 

(nt) Q = [1.203973  - ( 0 . 60xl06h) ( 2 . 0xl0_6f /h) ]/ ( 2 . 0xl0-6f /h) 

= 1987  unit-hours. 

The  POC  curve  for  this  plan  is  plotted  in  Figure  1. 

Step  7 : Let  us  first  examine  the  sensitivity  of  the  BAZE 

plan  to  changes  in  the  posterior  assurance  y and 
criterion  failure-rate  A^.  We  express  the  varying 
criterion  failure-rate  A^  as  a function  of  k 
according  to  A^  = kAg  = k(2.0xl0  ^f/h) . Table  I 
gives  the  resultant  test  plan  as  a function  of  a 
selected  grid  of  values  of  k and  y.  The  plan  (nt)g  = 
1987  unit-hours  is  indicated  in  table  for  k = 1.0 
and  y = 0.70.  It  is  clearly  observed  that  the 
required  unit-hours  of  test  increase  as  y increases 
and  k decreases.  It  is  observed  that  the  optimal 
test  plan  is  quite  sensitive  to  increasing  y and 
also  somewhat  sensitive  to  decreasing  k,  for  the 
range  of  k indicated.  The  "zeros"  indicate  those 
situations  in  which  the  prior  distribution  is 
sufficient  to  guarantee  that  the  risk  is  at  or 
below  the  specified  level,  without  the  need  for 
additional  testing. 
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Let  us  now  determine  the  sensitivity  of  the  BA2E 
plan  to  changes  in  the  prior  assurance.  Suppose 
that,  regarding  the  prior  distribution,  each  tail 
area  is  actually  2.5%,  rather  than  5%,  as  we  have 
assumed.  Applying  Steps  1-4  for  pQ  = 0.95  yields 
a = 1.36  and  b = 929.70xl03h.  Table  II  gives  the 
resultant  BAZE  test  plan  for  the  same  set  of  values 
of  y and  k used  in  Table  I.  Table  III  gives  the 
percentage  change  in  the  BAZE  plans  of  Table  II 
relative  to  the  "nominal"  plans  given' in  Table  I. 

It  is  observed  from  Table  III  that  the  BAZE  plans 
are  somewhat  sensitive  to  a 100%  error  (assumed  5% 
when,  in  fact,  2.5%)  in  the  tail  areas  of  the  prior 
distribution.  Now,  if  a = 1.36  and  b = 929.70xl0^h, 
the  "actual"  posterior  assurance  of  the  plan  (nt)g  = 
1987  unit-hours  is  easily  computed  from  (8)  to  be 
75%,  rather  than  70%  as  required.  Thus,  we  unknow- 
ingly would  have  more  posterior  assurance  than 
required.  Similarly,  for  any  other  BAZE  test 
plan  in  Table  I,  we  could  compute  the  "true"  posterior 
assurance  relative  to  this  "error"  in  fitting  the 
prior  distribution.  Similarly,  Tables  IV  and  V 
consider  the  case  where  the  actual  tail  areas  of  the 
gamma  prior  are  10%,  instead  of  5%,  as  assumed.  From 
Table  V,  it  is  observed  that  the  BAZE  plans  are  fairly 
insensitive  to  this  100%  error  in  the  prior  tail- area 


TABLE  IV.  Required  Unit-Hours  of  Testing  for  Selected  Values  of  k and 
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probabilities.  If  a = 0.67  and  b - 330.68xl03h,  the 
actual  posterior  assurance  of  the  plan  (nt) ^ = 1987 
unit-hours  is  66%,  as  compared  to  the  desired  value 
of  70%.  In  this  case,  we  have  less  posterior 
assurance  than  required. 

Step  8:  The  unconditional  probability  of  passing  the  test 

(nt)Q  = 1987  unit-hours  is  calculated  from  (10)  as 

P (Passing  the  Test)  = [ ( 0 . 60xl06 ) / ( 0 . 60xl06  + 1987)] 

= 0.9967. 

Table  VI  gives  the  unconditional  probability  of  pass- 
ing the  corresponding  test  given  in  Table  I.  In 
practice,  tables  such  as  Table  VI  are  useful  to  the 
producer  in  selecting  the  final  test  plan.  From  the 
results  of  Step  7 and  this  step,  the  final  plan  to 
be  used  is  (nt)^  = 1987  unit-hours  of  test. 

Step  9:  A single  component  may  be  tested  for  1987  hours; 

five  components  for  397.4  hours;  ten  components 
for  198.7  hours;  etc.  If  no  failures  occur,  it 
may  be  claimed  that  a failure-rate  of  2.0xl0~6f/h 
or  less  has  been  demonstrated  with  70%  assurance. 

At  this  point  the  question  can  be  raised;  namely,  how  many 
unit-hours  of  testing  are  required  when  using  an  alternate 
classical  (standard  non-Bayesian)  procedure?  By  judicious 
choice  of  producer  and/or  consumer  risks,  it  is  possible  to  "show" 


that  classical  test  plans  result  in  either  larger  or  smaller 
total  unit-hours  of  testing.  Thus,  a person  advocating  a purely 
classical  approach  could  "show"  that  his  procedure  results  in 
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less  testing.  On  the  other  hand,  a "Bayesian"  can  "show"  the 
opposite  to  be  the  case.  The  correct  answer  lies  in  recognizing 
that  the  two  procedures  cannot  be  directly  compared,  due  to  the 
basic  underlying  philosophical  difference  regarding  the  inter- 
pretation of  the  failure-rate  X.  In  the  classical  approach,  the 
failure-rate  is  a non-random  (unknown)  constant,  while  in  the 
Bayesian  approach  it  is  an  (unknown)  value  of  a random  variable. 
The  proper  procedure  to  be  used  should  be  based  on  this  fact 
alone,  and  not  on  the  basis  of  which  procedure  yields  the  smallest 
amount  of  testing.  In  some  cases,  the  classical  procedure  will 
require  less  testing;  in  others,  the  BAZE  procedure  will  require 
less.  Since  nuclear  reactor  components  exhibit  failure-rates 
which  appear  to  be  random,  a Bayesian  procedure,  such  as  described 
here,  should  be  used. 

VI.  PRIOR  DISTRIBUTIONS  AND  ANALYSIS  FOR  SELECTED  COMPONENTS 

OF  NUCLEAR  REACTOR  SAFETY  SYSTEMS 

Table  III  2-1  in  Appendix  III  of  the  WASH  1400  Reactor  Safety 
Study  (1975)  contains  assessed  estimates  of  the  failure-rates  of 
selected  components  of  PWR  and  BWR  safety  systems.  Upper  and 
lower  bounds,  as  well  as  the  median,  were  computed  from  best 
available  data  for  these  failure-rates.  The  approach  used  there 
was  also  a Bayesian  one  and  the  random  component  failure-rates 
were  assumed  to  follow  a log-normal  prior  distribution.  It  is 
remarked  here  that  the  log-normal  distribution  is  also  a two- 
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per  demand  and  per  hour  basis,  depending  upon  the  failure  mode, 
only  the  time-dependent  failure  modes  whose  failure-rates  are 
given  on  a per  hour  basis  will  be  considered  here.  The  upper 
and  lower  bounds  were  computed  to  be  5%  bounds,  respectively. 

That  is,  the  probability  that  the  given  interval  contains  the 
failure-rate  for  any  particular  system  is  90%.  Thus  p^  = 0.90. 
These  bounds  are  reproduced  in  Table  VII  for  these  same  components 
and  corresponding  failure  modes.  If  available,  failure-rates 
corresponding  to  US  nuclear  operational  experience  are  also  given. 

Gamma  prior  distributions  have  been  fitted  to  these  data  by 
means  of  the  procedure  outlined  in  Section  III.  These  prior  dis- 
tributions may  be  used  in  future  test  programs,  such  as  those 
required  in  demonstrating  the  safety  and  reliability  of  planned 
LMFBR  systems.  The  corresponding  gamma  parameters  a and  b,  as 
well  as  the  prior  mean  and  variance,  are  given  in  Table  VII.  It 
is  observed  that  there  are  only  three  different  "shapes"  of  gamma 
priors  present;  namely,  a = 0.50,  0.84,  and  2.45.  This  is  due  to 
the  fact  that  the  ratio  of  UL  to  LL  takes  on  only  three  different 
order  of  magnitude  values;  namely,  3,  2,  or  1.  The  priors 
corresponding  to  a = 0.50  and  0.84  are  quite  diffuse;  that  is, 
the  prior-variance  is  quite  large  for  these  priors. 

Table  VII  also  presents  certain  probabilities  of  interest 
associated  with  the  corresponding  fitted  gamma  prior  distributions. 
The  quantities  P^,  P2,  P^,  P^,  and  P^  in  Table  VII  represent  the 
probabilities  defined  by 
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5 Prob.  (A  < Prior  Mean), 

?2  = Prob.  (A  < Median  of  a Log-Normal  Prior), 

P^  = Proo.  (A  < US  Nuclear  Operating  Experience 
Value) , 

P4  = Prob.  (A  < Lower  5%  Limit)  , 

P^  = Prob.  (A  > Upper  5%  Limit) • 

These  probabilities  were  calculated  by  means  of  the  incomplete 
gamma  function  code,  INCGAM,  mentioned  in  Section  III.  Since  the 
values  of  P2  in  Table  VII  are  all  less  than  0.5,  the  gamma  prior 
distributions  tend  to  favor  somewhat  larger  failure-rates  than 
a log-normal  distribution  with  the  same  5%  tail- area  probabili- 
ties. The  gamma  approach  is  slightly  conservative.  Further,  the 
values  of  P4  and  P,-  are  included  as  a check  on  the  accuracy  of 
the  procedure  used  to  fit  the  gamma  priors.  The  P^  values  are 
all  5%,  to  three  decimal  places,  while  the  values  of  P^  deviate 
somewhat  from  the  desired  value  of  5%  for  a = 0.50  and  2.45.  This 
departure  is  due  to  errors  in  reading  and  interpolating  earlier, 
and  slightly  less  accurate,  versions  of  Figure  A1  and  Table  A2 , 
respectively . 

VII.  SUGGESTED  BAZE  TEST  PLANS  FOR  SELECTED  COMPONENTS  OF  NUCLEAR 
REACTOR  SAFETY  SYSTEMS 

The  purpose  of  this  section  is  to  propose  several  possible 
BAZE  test  plans  for  a few  typical  components  selected  from 
Table  VII.  It  is  emphasized  that  the  plans  given  here  are  not 
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Short  Across 
MO/NC  Contact 


(TEST)  Rupture 


Combination) 


to  be  considered  an  exhaustive  set  or  the  sole  source  for  the 
final  selection.  The  final  selection  rests  with  the  producer 
in  conjunction  with  the  consenting  agreement  of  the  consumer. 

It  may  be  necessary  or  desirable  to  alter  the  prior  distributions 
given  in  Table  VII  for  numerous  possible  reasons.  Such  altera- 
tions will  result  in  different  test  plans  from  those  proposed 
here. 

Before  such  BAZE  test  plans  can  be  used  in  practice,  several 
things  must  be  considered.  First,  the  life  cycle  of  each  compo- 
nent must  be  identified.  This  includes  establishing  the  environ- 
ments and  operating  conditions  to  which  each  component  is  exposed 
during  the  various  phases  of  its  existence.  This  also  includes 
defining  the  duty  cycle  for  the  operational  phases.  Such  practical 
details  surrounding  the  implementation  of  these  BAZE  plans  will 
not  be  addressed  here.  This  would  ordinarily  be  considered  in 
the  contractor's  reliability  program  plan. 

We  shall  restrict  our  consideration  to  three  "typical" 
components  given  in  Table  VII;  namely,  electric  clutchs  (a  = 0.84), 
electric  motors  (a  = 2.45),  and  the  category  of  flanges,  closures, 
and  elbows  (a  = 0.50^ . The  specified  failure-rate  X^  will  be 
taken  to  be  the  corresponding  value  of  X according  to  the  US 
nuclear  operating  experience  given  in  Table  VII.  If  such  values 
are  not  available,  as  indicated  in  Table  VII,  then  Xq  will  be 
taken  to  be  the  assessed  median  prior  failure-rate.  Once  a 
reliability,  risk,  or  safety  analysis  program  is  under  way,  the 
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value  of  Aq  is  a contract  specification  which  may  differ  from  the 
values  considered  here. 

First,  consider  the  electric  clutch  in  Table  VII.  The  BAZE 
procedure  described  in  Section  IV  was  applied  to  this  component 
in  conjunction  with  its  fitted  prior  given  in  Table  VII.  Table 
VIII  presents  the  resulting  required  unit-hours  of  testing  for 
selected  values  of  the  discrimination  ratio  (k)  and  the  posterior 
assurance  (y) • For  example,  for  k = 5 and  y = 0.90,  it  would 
require  (nt)  = [2.018547  - ( 26 . 72xl04 ) ( 5)  ( 1 . 0xl0~6 ) ]/[  ( 5 ) ( 1 . OxlO-6 ) ] 
136,510  unit-hours  of  electric  clutch  testing  without  a failure 
in  order  to  demonstrate  a failure-rate  of  5.0x10  ®f/h  with  90% 
assurance.  Table  IX  gives  the  unconditional  probabilities  of 
passing  the  corresponding  tests  given  in  Table  VIII.  The  uncon- 
ditional probability  of  passing  the  above  test  is  70.7%. 

Similarly,  Tables  X and  XI  give  the  required  BAZE  unit-hours 
of  testing  and  unconditional  probabilities  of  passing  the  test 
for  the  case  of  an  electric  motor  in  a normal  (nonextreme)  en- 
vironment. Likewise,  Tables  XII  and  XIII  consider  the  category 
of  flanges,  closures,  and  elbows,  with  respect  to  leak  or  rupture 
failure  modes. 

Now  let  us  consider  the  sensitivity  of  the  BAZE  procedure 
for  each  of  these  three  component  categories  to  errors  in  the 
tail- area  probabilities.  The  procedure  follows  the  example  in 
Section  V.  Table  XIV  gives  the  fitted  gamma  prior  distributions 
for  these  three  components  for  actual  prior  probabilities  of  being 
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TABLE  VIII.  Required  BAZE  Unit-Hours  of  Testing  for  the  Electric 
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TABLE  X.  Required  BAZE  Unit-Hours  of  Testing  for  the 


TABLE  XII.  Required  BAZE  Unit-Hours  of  Testing  for  Flanges, 
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below  and  above  LL  and  UL,  respectively,  the  tail-area  proba- 
bilities, of  2.5%,  5%,  and  10%.  Table  XV  gives  the  BAZE  test 
plans  for  the  electric  clutch  when  the  tail  area  probabilities 
are  2.5%  instead  of  5%.  Table  XVI  gives  the  percentage  change 
in  the  plans  of  Table  XV  relative  to  the  "nominal"  plans  of 
Table  VIII.  Similarly,  Tables  XVII  and  XVIII  correspond  to 
Tables  XV  and  XVI  except  that  the  tail-area  probabilities  are 
10%.  It  is  observed  from  Tables  XVI  and  XVIII  that  the  BAZE 
plans  for  the  electric  clutch  are  fairly  sensitive  to  the  tail- 
area  probabilities  in  the  gamma  prior  distribution.  Similarly, 
Tables  XIX-XXII  and  Tables  XXIII-XXVI  correspond  to  Tables 
XV-XVI I I for  the  case  of  the  electric  motor  and  flanges,  closures, 
and  elbows,  respectively.  Again,  the  BAZE  plans  for  these  com- 
ponents are  observed  to  be  fairly  sensitive  to  the  tail-area 
probabilities. 

Finally,  similar  BAZE  test  plans  could  be  developed  for  the 
remaining  components  in  Table  VII.  The  BAZE  procedure  appears 
to  be  a practical  means  of  conducting  Bayesian  reliability  de- 
monstration testing.  The  main  advantage  over  other  Bayesian 
procedures  is  its  ease  of  application.  A few  simple  tables  and 
figures,  which  are  included,  and  a pocket  calculator  are  all 
that  is  required.  The  practical  utility  of  the  procedure 
appears  to  be  evident. 
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TABLE  XIV.  Gamma  Prior  Distributions  for  Selected  Tail-Area  Probabilities  for  the 
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Relative  to  Table  VIII 
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TAME  XVIII.  Percentage  Change  In  the  HAZE  Plans  of  Table  XVII  Relative  to  Table  VIII 


TABI.E  XXII.  Percentage  Change  in  the  BAZE  Plana  of  Table  XXI  Relative 
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APPENDIX 

Figure  Al . 
Figure  A2. 

Table  Al . 
Table  A2. 


. Some  Useful  Figures  and  Tables  for  the  BAZE  Procedure 

Gamma  Shape  Parameter  a as  a Function  of  log^g  (UL/LL) 

Gamma  Reference  Scale  Parameter  bg  As  a Function  of 
the  Shape  Parameter  a 

Values  of  0 for  Selected  Values  of  Prior  Shape 
Y 

Parameter  (a)  and  Posterior  Assurance  (y) 

Gamma  Reference  Scale  Values  bg  for  Selected  Values 
of  Shape  Parameter  (a)  and  Prior  Assurance  (pg) 
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Table  Al.  Values  of  e y for  Selected  Values  of  Prior  Shape  Parameter  (a) 

and  Posterior  Assurance  00 
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TABLE  A2 . 


Gamma  Reference  Scale  Values  bg  for  Selected  Values 
of  Shape  Parameter  (a)  and  Prior  Assurance  (pg) 


a/p0 

0.95 

0.90 

0.80 

0.01 

3-5227-155 

4.4655-125 

5.6607E-95 

0.02 

4.5019E-75 

5.0687E-60 

5.7063E-45 

0.03 

2.2797E-48 

2.4673E-33 

2.6702E-23 

0.04 

5. 1503E-35 

1 .7282E-27 

5.7988E-20 

0.05 

5 . 3157E-27 

5.5739E-21 

5 . 8446E- 1 5 

0.06 

1 .1726E-21 

1 .2199E-16 

1 . 269  IE-1 1 

0.07 

7.7079E-18 

1 . 5394  E— 1 3 

3.07U6E-09 

0.08 

5.6375E-15 

3-2656E-1 1 

1 .891 6E-07 

0.09 

9.5333E-13 

2. 1098E-09 

4.6669E-06 

0.10 

5.7917E-11 

5.9307E-08 

6.0730E-05 

0.11 

1 .669CE-09 

9. 1013E-07 

4.9629E-04 

0.12 

2.7502E-08 

8 .8704E-06 

2.861  IE-03 

0.13 

2.9479E-07 

6.0971E-05 

1 .261  IE-02 

0. 14 

2.2538E-06 

3.1852E-04 

4.501 4E-02 

0.15 

1 .3150E-05 

1 . 3360E-03 

1 .3573E-01 

0.16 

6. 1598E-05 

4.6882E-03 

3 . 5 682E-0 1 

0.17 

2.4079E-04 

1 .4204E-02 

8.3786E-01 

0.18 

8 .0955E-04 

3 . 8074E-02 

1 .7907E+00 

0.19 

2 . 3972E-03 

9 . 2062E-02 

3.5355E+00 

0.20 

6 . 3725E-03 

2.0392E-01 

6 .5255E+00 

0.21 

1 .5443E-02 

4.  1900E-01 

1 .13635+01 

0.22 

3.4551E-02 

3 . 0632E-01 

1 .3841E+01 

0.23 

7.21 12E-02 

1 .4684E+00 

2.99005+01 

0.24 

1 .4163E-01 

2.5435E+C0 

4.56825+01 

0.25 

2.6366E-01 

4 . 2 1 86E+00 

6.75015+01 

0.26 

4.6315E-01 

6 . 7327E+00 

9.6S35E+01 

0.27 

7 . 9 69  7 E— 0 1 

1 .03S4E+01 

1 . 3531E+02 

0.28 

1 .3067E+00 

1 .5534E+01 

1 .8469E+02 

0.29 

2 .0715E+00 

2.261 1E+01 

2.4635E+02 

0.30 

3. 1857E*00 

3.21 105+01 

3.2372E+02 

0.31 

4.7668E+00 

4 . 4597E+01 

4. 173UE+O2 

0.32 

6.9577E+00 

6 . 0702E+0 1 

5.29765+02 

0.33 

9 92835+00 

8.1121 E+0 1 

6 . 6303E+02 

0.34 

1.38805+01 

1 .C661E+02 

3.1924E+C2 

0.35 

1 .9041E+01 

1 . 379SE+02 

1 .C0C4E+03 

0.36 

2.5675E+01 

1.76105+02 

1.20365+03 

0.37 

3.4Q74E+01 

2 . 2 1365+02 

1 . 4U57T+C3 

0.38 

4 . 4566E+0 1 

2 . 7622E+G2 

1.71 35E*Q3 

0.39 

5 . 7506E+01 

3 .40155+02 

2.01 4CE+C  3 

0.40 

7 . 3283E+0 1 

4. 1465E+02 

2 .3489E+03 

0.41 

9 .2314E+01 

5.CC7ME+02 

2 . 7 1 97£+0  3 

0.42 

1 . 1504E+02 

5.9945E*02 

3.12305+03 

0.43 

1.419UE+02 

7 . 1 1 80E+02 

3.57525+03 

0.44 

1 .7351E+02 

3.3884E*02 

4.0627E+03 

0.45 

2 . 1 026E+02 

9.31605*02 

4.591 7E+03 

0.46 

2.52725+02 

1 . 141 1E*03 

5.  1633E+03 

0.47 

3 . 01 46E+02 

1 . 3 1 S3E*03 

5.77865*03 

0.48 

3 • 5703E+02 

1.51425*03 

6.43355+03 

0.49 

4.20025+02 

1.72985*03 

7.14  38  5+0  3 

0.50 

4.9103E+02 

1 .96615*03 

7.8954E+03 

0.51 

5 .7065E+02 

2.2 2385*03 

8.6938E+03 

0.52 

6 . 5947E+02 

2 . S039E+03 

9 .5397E+03 

0.53 

7.531 0E+02 

2 . 30732*03 

1 .043UE+CU 

0 . 5« 

3 . 67 1 5E+C2 

3.13475+03 

1 . 13755+04 

0.55 

9 • 8723E+02 

3.4369E+03 

1 .23675+04 

0.56 

1.1 1 89E*03 

3.8647E+03 

1 . 34Q7E+Q4 

0.57 

1.26235+03 

u. 26835+03 

1 . U496E+04 

0.58 

1.41965+03 

4 . 6997E+03 

1 . 56355+OU 

0.59 

1.5397E*03 

5.15325*03 

1 .6823E+04 

0.60 

1 .7727E*03 

5 . 64U8E+03 

1 .306CE+04 

Table  A2  (Continued) 


«/.P5 

0.95 

0.90 

0.61 

1.9723E+03 

6 . 1 602E+03 

0.62 

2. 1859E+03 

6 . 7047E+03 

0.63 

2 . 4 1 52E+03 

7.2789E+03 

0.64 

2 .66058+03 

7 . 8833E+03 

0.65 

2.92258+03 

8.51 82E+03 

0.66 

3 . 2016E+03 

9 . 1840E+03 

O.67 

3.4982E+03 

9.381 1E+03 

0.68 

3.81 30E+03 

1 .G610E+04 

0.69 

4 . 1 462E+03 

1 . 1370E+04 

0.70 

4.4983E+03 

1 . 2 1 63E+04 

0.71 

4.8698E+03 

1 .2983E+04 

0.72' 

5.261 OE^O  3 

1 .3846E+04 

0.73 

5.6722E+03 

1 .4737E+04 

0.74 

6. 1040E+03 

1 .5660E+04 

0.75 

6 . 5565E+03 

1 .6616E+04 

0.76 

7 . 0302E+03 

1 .7606E+04 

0.77. 

7 .5253E+03 

1 .8629E+04 

0.78 

8.0421E+03 

1.96358+04 

0.79 

8.58098+03 

2 . 0775E+04 

0.80 

9.1420E+03 

2. 1898E+04 

0.81 

9 .7256E+03 

2. 3C54E+04 

0.82 

1 .0332E+04 

2.4244E+04 

0.83 

1 .0961 E+04 

2.5467E+04 

0.84 

1.1614E+04 

2.6723E+04 

0.85 

1 .2289E+04 

2 . 801 3E+G4 

0.36 

1 .2989E+04 

2.9336E+04 

O.87 

1 .3712E+04 

3 . 0692E+04 

0.88 

1 .4459E+04 

3 . 2082E+04 

0.89 

1 . 5229E+04 

3 . 3504E+04 

0.90 

1 .6024E+04 

3 • 4959E+04 

0.91 

1.6843E+04 

3. 6448E+04 

0.92 

1 .7687E+04 

3.7968E+04 

0.93 

1 .8555E+04 

3 . 9522E+04 

0.94 

1 .9447E+C4 

4.11 07E+04 

0.95 

2.0364E+04 

4 . 2725E+04 

0.96 

2.1 305E+04 

4 . 4 375E+04 

0.97 

2 . 227 1 E+04 

4 . 6057E+04 

0.98 

2 . 3262E+04 

4.7771E+04 

0.99 

2.4278E+04 

4.951  6E+C4 

1.00 

2.5318E+04 

5 . 12938+04 

1.05 

3 . 089 1 E*C4 

6.0643E+04 

1.10 

3.7082E+04 

7 . 0752E+04 

1.15 

4 . 3383E*04 

3. 1594E+04 

1.20 

5.1301 E+04 

9 . 31458+04 

1.25 

5.9311 E+04 

1 .0538E+05 

1.30 

6.7903E+04 

1 . l827E+05 

1.35 

7.7079E+04 

1.31798+05 

1.40 

8.6310E+04 

1.45938*05 

1.45 

9.7087E*04 

1 .6064E+05 

1.50 

1 .0790E+05 

1 .7592E+05 

1.55 

1 .19238+05 

1 .9174E+05 

1.60 

1.3106E+05 

2 . Q803E+05 

1.65 

1 .4338E+05 

2 . 2492E+05 

1 .70 

1 .5618E-05 

2 . 42258*05 

1.75 

1 .6944E+05 

2 . 60038*05 

1 .80 

1 .8314E+05 

2 . 7827E+G5 

1.85 

1 .9729E-.05 

2.96938*05 

1.90 

2. 1186E+05 

3. 1601E+05 

1.95 

2.2633E+05 

3.35498*05 

2.00 

2.4221 E+05 

3.5536E*05 

0.80 


1 .9348E+04 
2.0685E+G4 
2 . 2071 E*04 
2 . 3506E+04 
2.4991E+04 
2.6525E+04 
2 . 8 1 08E+04 
2 . 9739E*04 
3. 1418E+04 
3 . 3 1 45E+04 
3.4921 E+04 
3 . 6743E+04 
3.8612E+04 
4.0523E+04 
4.2491 E+04 
4.4499E+04 
4 . 6552E+04 

4.8651 E+04 
5 • 0794E+04 
5.29S2E+04 
5.521 3E+04 
5.7488E+04 
5 . 93C6E+04 
6.2166E+04 
6.4569E+04 
6.701 3E+04 
6.9493E+04 
7 . 2024E+04 
7.4590E+04 
7.7197E*04 
7. 9842E+04 
8 . 2527E+04 
8.5251 E+04 
8.801 2E+Q4 
9.081 2E+04 
9.3649E+04 
9 . 6522E+G4 
9 . 9432E+04 

1 .0238E+05 
1 .0536E+05 
1 . 20798+05 
1 .3705E-05 
1 .5410E+05 
1 .719nE*05 
1 .9039E+05 
2.0956E+05 
2.2 935E+05 
2.497UE+05 
2.7069E+05 
2.921 9 Ei-05 
3. 1419E-05 
3. 3669E-05 
3 • 5965E+05 
3.8306E+C5 
4 . 0689E+05 
4.31 13E*05 
4.5575E+05 
4.8075E+C5 
5.061 1E*05 
5.3181 E*05 
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0.95 


0.90 


0.80 
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Table  A2  (Continued) 


a/p0 


2.05 

2.5797E+05 

3.7560E+05 

5.5784E+05 

2.10 

2.7411 E+05 

3-9621E+05 

5 . 34  igE+05 

2.15 

2 .9062E+05 

4. 1716E+05 

6 . 1 085E+Q5 

2.20 

3 . 0748E+05 

4.3845E+05 

6 . 3780E+05 

2.25 

3 .2469E+05 

4.6007E+05 

6.6504E+05 

2.30 

3.4224E+05 

4.8201 E+05 

6 . 9255E+05 

2.35 

3.6011 E+05 

5 . 0425E+05 

7 . 2032E+05 

2.40 

3.7830E+05 

5 . 2679E+Q5 

7.4835E+05 

2.45 

3 - 9680E+05 

5 . 4963E+05 

7.7663E+05 

2.50 

4.1561 E+05 

5.727  4E+05 

8 . 051 5E+05 

2.55 

4.3470E+05 

5 .9612E+05 

8.3391E+05 

2.60 

4 . 5409E+05 

6 . 1977E+05 

8 . 6288E+05 

2.65 

4.7375E+05 

6 .4363E+05 

8.9208E+05 

2.70 

4.9369E+05 

6 . 6784E+05 

9 . 21 49E+05 

2.75 

5 . 1 389E+05 

6.9224E+05 

9.51 10E+05 

2.80 

5 . 3435E+05 

7. 1688E+05 

9 . 8092E+05 

2.85 

5.5507E-05 

7 . 4175E+05 

1 .01C9E+06 

2.90 

5.7603E+05 

7 . 6635E+C5 

1 .041 1E+06 

2.95 

5.9723E+05 

7 . 92 1 6E+05 

1 .0715E+06 

3.00 

6. 1867E+05 

8 . 1769E+05 

1 . 1 021 E+0  6 

3-05 

6.4034E+05 

8.4343E+05 

1 . 1328E+06 

3.10 

6 . 6224E+05 

8 . 6937E+05 

1 . 1637E+06 

3-15 

6.8435E+05 

8.9551 E+05 

1 . 1948E-06 

3-20 

7.0668E+05 

9 . 2184E+05 

1 .2260E+06 

3-25 

7 . 2922E+05 

9.4837E+G5 

1 . 257  4E+C6 

3-30 

7.5197E+05 

9 . 7507E+05 

1 .2889E+C6 

3-35 

7.7492E+05 

1 .0020E+06 

1 .3206E+06 

3-40 

7.9807E+05 

1 .0290E+06 

1 .3525E+06 

3-45 

8.2141 E+05 

1 .0563E+06 

1 .3844E+06 

3-50 

8 . 4493E+05 

1 .0837E+06 

1 .4166E+06 

3-55 

8 . 6865E+05 

1.1112E+C6 

1 .4438E+06 

3.60 

8.9255E+05 

1 . 1390E+06 

1 .4S12E+06 

3-65 

9 . 1662E+05 

1 . 1669E+06 

1 .5137E+06 

3.70 

9 . 4087E+05 

1 . 1949E+06 

1 .5464E+06 

3-75 

9 . 6529E+05 

1 .2231E+06 

1 .5791E+06 

3-80 

9.3988E+05 

1 .2515E+06 

1 .6120E+06 

3-85 

1 .0146E+06 

1 .2800E+06 

1 .6450E+06 

3-90 

1 .0396E+06 

1 .3086E+06 

1 .6782E+06 

3-95 

1 .0646E+06 

1 .3374E+06 

1 .711 4E+0  6 

4.00 

1 .0899E+06 

1 .3663E+06 

1 .7448E+06 

4.05 

1.1153E+06 

1 .3954E+06 

1 .7782E+06 

4.10 

1 . 14C8E+06 

1 .4246E+C6 

1 .31 18E+06 

4.15 

1 . 1665E+06 

1 .4539E+06 

1 .8455E+06 

4.20 

1 . 1923E+06 

1 .4333E+06 

1 .8793E+06 

4.25 

1 .2183E+06 

1 .5129E+06 

1 .9132E+C6 

4.30 

1 .2444E+06 

1 .5426E+06 

1 .9472E+06 

4.35 

1 .2706E+06 

1 .5724E+06 

1 .9812E+06 

4.40 

1 .2970E+06 

1 .6023E+06 

2 . 0 1 54E+06 

4.45 

1 • 3235E+06 

1 .6324E+06 

2 . 0497E+06 

4.50 

1 . 3502E+06 

1 . 66262^06 

2.0841 E+06 

4.55 

1 . 3770E+06 

1 .6928E+06 

2 . 1 135E+06 

4.60 

1.4039E+06 

1 .7232E+06 

2.1531 E+06 

4.65 

1 .4309E+06 

1 .7537E+06 

2. 1877E+06 

4.70 

1 . 458 1 E+06 

1 .7343e*06 

2.222 5E+06 

4.75 

1.4853E+06 

1 .8 1 50E+06 

2.2573E-06 

4.30 

1 .5127E-06 

1.3459E-06 

2.2922E+06 

4.35 

1 . 5403E+06 

1 .8768E*06 

2 . 3272E+06 

4.90 

1 .5679E+06 

1 .9079E*C6 

2.’622E+06 

«.95 

1 . 5956E+06 

1 .9389E+06 

2.3974S*06 

5.00 

1 .6235E+06 

1 .9701S+C5 

2 . 4326E*06 

Table  A2  (Continued) 


*/po 

0.95 

0.90 

0.80 

5.05 

1.6515E-f06 

2.0015E+06 

2.4679E+06 

5.10 

1 .6795E+06 

2 . 0329E+06 

2.5033E+06 

5.15 

1 .7077E+06 

2.0644E+06 

2.5387E+06 

5.20 

1 .7360E+06 

2.0960E+06 

2 . 5742E+06 

5.25 

1 .7644E+06 

2. 1277E+06 

2.6098E+06 

5.30 

1.7929E+06 

2. 1594E+06 

2.6455E+06 

5.35 

1.8215E+06 

2.191 3E+06 

2.6813E+06 

5.40 

1.8502E+06 

2 . 2232E+06 

2.7171 E+06 

5.45 

1 .8790E+06 

2 . 2553E+06 

2.7529E+06 

5.50 

1.9079E+06 

2.2874E+06 

2.7839E+06 

5.55 

1.9369E+06 

2. 3196E+06 

2.8249E+06 

5.60 

1.9659E+06 

2.3519E+06 

2.861 0E+06 

5.65 

1.9951E+06 

2 . 3843E+06 

2.8971E+06 

5.70 

2.0244E+06 

2 . 4 1 67E+06 

2 . 9333E+06 

5.75 

2.0533E+06 

2. 4492E+06 

2 . 9696E+06 

5.80 

2.0832E+06 

2. 4818E+06 

3.0060E+06 

5.85 

2.11 28E+06 

2 . 5 1 45E+06 

3 . 0423E+06 

5.90 

2. 1424E+06 

2 . 5473E+06 

3 . 0788E+06 

5.95 

2.1721E+06 

2.580 1E+06 

3. 1153E+06 

6.00 

2.2019  E+06 

2.6130E+Q6 

3. 1519E+06 

6.05 

2 . 2318E+06 

2 . 6460E+06 

3. 1885E+06 

6.10 

2 . 2617E+06 

2 . 6790E+06 

3 . 2252E+06 

6.15 

2.2918E+06 

2.7122E+06 

3 . 2620E+06 

6.20 

2 . 3219E+06 

2.7454E+06 

3.2988E+06 

6.25 

2.3521E+06 

2.7786E+06 

3-3356E+06 

6.30 

2.3824E^06 

2.81 19E+06 

3-372-6  E+06 

6.35 

2 . 41 28E+06 

2.8453E+06 

3 . 4095E+Q6 

6.40 

2.4433E+06 

2.8788E>06 

3 . 4465E+06 

6.45 

2.4738E+06 

2 .9  123E+06 

3 . 4836E+06 

6.50 

2 . 5044E+06 

2.9459E+06 

3.5208E+06 

6.55 

2.5351E+06 

2.9796E+06 

3.5579E-h06 

6.60 

2 . 5658E+06 

3 . 01 33E+06 

3 . 5952E+06 

6 . 65 

2. 5966E+06 

3 . 0471 E+06 

3 . 6324E+06 

6.70 

2.6275E+06 

3.0810E+06 

3 . 6698E>06 

6.75 

2. 6585E+06 

3. 1 149E+06 

3 .7072E+06 

6.80 

2 . 6895E+-06 

3 . 1 488E+06 

3.7446E+06 

6.85 

2.7206E+06 

3 • 1829E+06 

3.7821E+06 

6.90 

2 .7518E+06 

3 .2170E+06 

3.3196E+06 

6.95 

2.7831E+06 

3.251 1E+06 

3 . 8572E+06 

7.00 

2.81 44E+06 

3 . 2853E+06 

3.8948E+06 

7.05 

2 . 8457E+Q6 

3.3196E+06 

3 .9324E+06 

7.10 

2 . 8772E+Q6 

3 • 3539E+06 

3.9701E+06 

7.15 

2.9087E+06 

3 . 3883E+06 

4.0079E+06 

7.20 

2.9403E>06 

3 . 4227E+06 

4.0457E+06 

7.25 

2.9719E>06 

3 . 4572E+06 

4.0835E+06 

7.30 

3.0036E+06 

3.4918E+06 

4 . 1 2 1 4E+06 

7.35 

3.0354E+06 

3.5264E+06 

4 . 1 593E+06 

7.40 

3 . 0672£fQ6 

3.5610E+06 

4 . 1973E+06 

7.45 

3.0991E+06 

3 .5957E+06 

4 . 2353E+06 

7.50 

3. 131 1E+06 

3.6305E-rOo 

4.2734E+06 

7.55 

3. 1631E>06 

3 .6653E+06 

4.31 15E+06 

7.60 

3 . 1952E+06 

3 . 7001 E+06 

4 . 3496E-*-06 

7.65 

3.2273E+06 

3.7351E+06 

4 . 3878E+06 

7.70 

3 • 2595E+06 

3 .7700E+06 

4 . 4260E+06 

7.75 

3 .29 17E+-06 

3.3050E+06 

4.4643E+06 

7.80 

3.3240E>06 

3.3401E+06 

4 . 5026S+06 

7.85 

3 • 3564E+06 

3 . 8752E*06 

4 . 5409E+06 

7.90 

3 . 3888E+06 

3.910uE+06 

4 . 5793E*Q6 

7.95 

3 . 421 3E+06 

3.9456E-06 

4.61 77E+06 

8.00 

3 . 4538E+-06 

3.9808E+06 

4 . 656 1 E*0  6 
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«/p0 

0.95 

0.90 

8.05 

3 . 4864E+06 

4.0161  E-^06 

8.10 

3.5191E+06 

4.0515E+06 

8.15 

3-5518E+06 

4 . 0869E+06 

8.20 

3 . 5845E+06 

4. 1223E+06 

8.25 

3.6173E+06 

4 . 1578E+06 

8.30 

3 . 6502E+06 

4 . 1933E+06 

8.35 

3.6831E+06 

4.2289E+06 

8.40 

3 . 71 60E+06 

4.2645E+06 

8.45 

3.7490E+06 

4 . 3002E+06 

8.50 

3.7821E+06 

4 . 3359E.+06 

8.55 

3.8152E+06 

4.3716E+06 

8.60 

3 . 8484E+06 

4.4074E+06 

8.65 

3.8816E+06 

4 . 4433E+06 

8.70 

3 .91 «8E+06 

4.4791E+06 

8.75 

3 .9481E+06 

4.5150E+06 

8.80 

3 .9815E+06 

4.551 0E+06 

8.85 

4.01  49E+-06 

4 .5870E+06 

8.90 

4 .0483E+06 

4.6230E+06 

8.95 

4.0818E+06 

4 . 6591 E+06 

9.00 

4 . 1154E+06 

4.6952E+06 

9-05 

4 . 1490E+06 

4.731 4E+06 

9.10 

4. 1826E+06 

4 . 7676E+06 

9-15 

4.2163E+06 

4.8038E+06 

9.20 

4 . 2500E+06 

4.8401 E+06 

9.25 

4.2838E+06 

4 . 8764E+06 

9.30 

4 . 3176E+06 

4.91 27E+06 

9-35 

4.351 4E+06 

4.9491E+06 

9.40 

4.3853E+06 

4.9856E+06 

9.45 

4 . 4193E+06 

5 . 0220E+06 

9-50 

4.4533E+06 

5.05S5E+06 

9.55 

4.4873E+06 

5.0950E+06 

9.60 

4.521 4E+06 

5 . 1316E+06 

9-65 

4 . 5555E+06 

5 . 1 682E+06 

9.70 

4 . 5896E+06 

5 .2048E+06 

9.75 

4 .6238E+06 

5.2415E+06 

9.80 

4 . 6580E+06 

5 . 2782E+-06 

9.85 

4 . 692  3E+06 

5.3150E+06 

9-90 

4 . 7266E+06 

5.3518E+06 

9-95 

4 .7610E+06 

5.3886E+06 

10.00 

4 . 7954E+O6 

5.4254E+06 

0.80 


4.6946E+06 
4.7331E+06 
4.7717E+06 
4.3103E+06 
4.8489E+06 
4.8876E+06 
4.9263E+06 
4.9650E+06 
5.0038E+06 
5 . 0426E+06 
5 . 08l 4E>06 
5 . 1203E+06 
5. 1592E+06 
5.1982E+06 
5 .2371E+06 
5 .2761E+06 
5 . 3152E+06 
5 . 3542E+06 
5 . 3933E+06 
5 . 4325E+06 
5.471 6E+06 
5.51 08E+-06 
5.5500E+06 
5 . 5893E+06 
5 . 6286E+06 
5 . 6679E-t-06 
5 .7072E+06 
5 .7466E+06 
5 . 7860E+06 
5 .8255E+06 
5.8649E+06 
5 . 9044E+06 
5.9439E+06 
5 . 9835E+06 
6.0230E+06 
6 . 0626E+06 
6 . 1023E+06 
6. 1419E+06 
6. 181 6E-*-06 
6 . 221 3E+06 


APPENDIX  B.  A Procedure  for  Selecting  a Gamma  Prior 
Distribution 


It  is  required  to  find  a and  b which  satisfies 


UL 


/ 

LL 


ba 
r (a) 


a-1  -bx. 
x e dx 


Po 


t 


(Bl) 


where 


LL  , a . . c 

r b a-1  -bx,  c 

L flaT  e dx  - / 


UL 


b a-1  -bx , . , . 

TUT  x e dx  = (p-Pq^2' 


Letting  y = x/LL  in  (Bl) , we  have 


fUL/LL  (bLL) a a-1  - (bLL) y 
\ r (a)  Y e 


dY  = P0* 


(B2) 


Since  b is  a scale  parameter,  set  bLL=l,  and  solve  (B2)  for 
the  shape  parameter  a.  Thus  a depends  only  upon  the  value  of  UL/LL, 
or  equivalently,  log(UL/LL).  Once  a has  been  numerically  determined, 
we  can  solve  (B2)  for  a temporary  value  of  b,  say  bQ,  corresponding 
to  a temporary  lower  limit  of,  say,  1.0x10  6f/h.  Since  b is  a 
scale?  parameter,  we  know  that 

bLL  = bQ ( 1 . OxlO-6 ) 

from  which 


b = bQ (1 . 0x10  6f/h) /LL. 
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